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Abstract 

The scheme for measurement of the state of a single spin (or a few spin system) based 
on the single-electron turnstile and injection of spin polarized electrons from magnetic metal 
contacts is proposed. Applications to the recent proposal concerning the spin gates based 
i^J ' on a silicon matrix (B.Kane, Nature, 393, 133 (1998)) are discussed. 

Q^ ■ After the discovery of efficient quantum algorithms [1] and a rigorous proof of the possibility 

of fault-tolerant quantum computing [2] , various realizations of quantum logical gates were pro- 
posed based on cold ions [3] , nuclear magnetic resonance [4] , optical schemes [5] , semiconductor 

I ^ ■ heterostructures [6], and Josephson effect [7]. Recently, the possibility of fabrication of quantum 

gates employing a silicon matrix with the dopant P 31 atoms was suggested [8] . The role of quan- 

£S) | turn bits is played by the nucleus and electron spins of the P 31 atom. One of the problems in the 

scheme [8] is the measurement of a single nucleus or electron (or both) spin. The papers [9,10] 
consider the indirect spin measurement employing the single-electron transistor. However, the 
proposed schemes do not allow to measure the state of a single spin and only make it possible 
to measure the different charge states of a system of nucleus and/or electron spins. 

The detection of a single spin in itself is not an exotic thing. The observation of Larmor 
precession of a single spin localized on the Si(lll)7x7 surface in ultra-high vacuum with the 
scanning tunnelling microscope (STM) was first reported by the IBM group (Demuth et al [11], 
see also [12]) ten years ago. There were also reports on the detection of the electron paramagnetic 
resonance signal in the STM current from a single spin in an organic molecule [13]. The STM 
with a magnetic tip was also demonstrated to be sensitive on the atomic level to the state of 
single spins on the surface of magnetic materials [14]. The quantity which is directly measured 
in STM is the tunnel current which depends on the lateral tip position relative to the sample 
on the atomic scale, and in the case of a magnetic tip the tunnel current has a spin-dependent 
^ ' component [12] 

y~J ■ ^t(x) oc p c {~x)pt m c (x) • m t , 

where p c ,t( x ) are the local densities of states on the tip apex and the sample surface, respectively, 
while m c (x) and mt are the local magnetizations at the sample surface at a point x and the tip 
apex, respectively. However, the steady state current measurements cannot be directly applied 
to the detection of the state of quantum gates. 

Realization of quantum gates requires the possibility of measuring the state of the system at 
an arbitrary moment of time. According to the general theory of quantum-mechanical measure- 
ments [15-17], the most complete description of any particular measuring procedure which can 
be applied to a quantum system is given by the so-called instrument. The instrument T{d\) is 
actually a mapping of the set of all quantum states of the system (density matrices) p s before 
the measurement to the system states (up to the normalization) just after the measurement 
p s = T(d\)p s which gave the result in the interval dX, the probability of obtaining a measure- 
ment result in the interval d\ being Trp s = Ti{T(dX)p s }. Any instrument can be represented 
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in the form 

T{d\) Ps = Tt a {(I s ® M A {d\))U(p s ® PaW' 1 }, 

i.e. any measurement can be realized by allowing the studied system to interact with a suitable 
auxiliary system (ancilla) prepared in a fixed initial state pa for some time (so that their joint 
evolution is described by certain unitary dynamics U) and a subsequent measurement generated 
by a suitable identity resolution M^(dA) performed on the ancilla. 

Proposed below is the method for detection of the state of a single spin (or a few spins, e.g. 
nucleus spin + electron spin) based on the "turnstile" concept [18,19]. Explicitly present in the 
proposed scheme are the preparation of the ancilla (pa) at an arbitrary moment of time, turning 
on the interaction between the ancilla and the studied system, their joint unitary evolution, 
turning off the interaction at an arbitrary moment of time, and detection of the state of the 
ancilla. Consider the following model system. Suppose that we have a system of spins in the 
subsurface region, e.g. an atom with the nucleus possessing non-zero spin and an electron 
localized on that atom (Fig. 1). Suppose also that a system of tunnel-coupled quantum dots is 
fabricated on the surface in such a way that the central dot is located just above the system of 
spins acting as the quantum bits (Figs. 1,2). Each dot has a single size-quantized level. The 
leftmost and the rightmost dots are tunnel-coupled to the magnetic metal electrodes (Figs. 1,2). 

The Hamiltonian describing the interaction among the quantum dots and between the quan- 
tum dots and metallic electrodes can be written as 

H = Yl £ ka a at<ra a k<ra + X^ £cC ^ CcCT + £ L c La C La + S R C R \ T C Ra )+ (1) 

k,<r,a=L,R & 

Y{TkLc\ a a kaL + T cL c+ a c La + T cR c+ a c Ra + h.c.) + Y(U L n La n L _ a + \J c n C(J n c ^ (J U r n Ra n R - (y ), 
kcr a 

where the first two terms describe the electron states in isolated electrodes and the dots, the 
third one represents the tunnel coupling between the dots and the electrodes, and the last term 
accounts for the Coulomb intradot repulsion (if it is important). We assume that the electrons 
in the electrodes are spin-polarized with the magnetization vectors (n^ and n R in the left and 
right electrodes, respectively) fixed by e.g. magnetic anisotropy. If the system is placed in an 
external magnetic field, the corresponding terms should be added to the Hamiltonian. The spin 
system (quantum gate) Hamiltonian, for example, for the case of nucleus spin + the electron 
spin localized on it can be written as 

H S = Y^ £ s ct a c sa + g s pBct a c sa/ a a(T ' ■ H + ginsl ■ H + g sI l ■ a aa < c+ a c S(7 > . (2) 

(j 

In the external magnetic field the contribution from the metal electrodes should also be taken 
into account. The Hamiltonian describing the interaction of the spins in the quantum gate and 
the electron localized in the central dot (see below) depends on the specific geometry of the 
considered structure. For example, if the wave functions of the central dot electron and the 
electron localized on the subsurface center overlap, then the Hamiltonian can be written as 

Hint = ^2(t sc C+ a C ca + h.C.) + g cI I ■ (Tva'C^Cca' + ^ Uscn C an s -a- (3) 

(j a 

If the overlap is negligible, then only the dipole-dipole interaction should be retained. 



The complete solution of the problem of finding the temporal evolution of the considered 
system is a difficult task. To go one step further, we shall assume that the characteristic times 
of different processes occurring in our systems form a certain hierarchy. Let T res be the typical 
time of electron tunneling from the central dot to the metal electrode when the energy levels 
in the adjacent quantum dots are tuned to the resonance (this time actually coincides with the 
time required for the electron to tunnel from the left (right) quantum dot to the metal electrode 
through a single barrier), T non be the characteristic tunneling time from the central dot to the 
metal electrodes when the levels in the dots are detuned from the resonance, and finally Td yn be 
the typical time of joint evolution caused by the interaction between the electron in the central 
dot and the spins in the gate. We shall assume that r res <C Td yn <C r non . Below we wish to take 
advantage of the well known point that for the tunneling through the two barriers (from the 
central dot to the metal electrodes), tuning of the levels into the resonance lifts the smallness 
associated with the additional barrier. The characteristic time are inversely proportional to the 
level width and depends on the level position; an estimate is given by (e.g., see Ref.[19]) 

"TT = ^) = r( ^ TL ~ C? <%+- 2 ' 70 = £ |TfeL|25(w " £kL) = |Tl|2 - 
r(w) [e c {Lo)-e L (uj)] 2 +^ Y 

Here 70 = \Tl\ 2 ~ |?cZ,| 2 = \T\ 2 is the bare tunnel transparency of the barrier between the dots 
and between the dots and the electrodes which can be assumed to be the same without any loss 
of generality. Under the resonance conditions (e c (uj r ) = £x(av)) ^/T res ~ \T\ 2 = 70- When the 
levels are detuned by the energy larger than the level width (A 3> 70), the characteristic time 
becomes l/r non « 70 (70 /A 2 ) <C l/r res . Accounting for the Coulomb repulsion does not change 
the situation qualitatively. 

Let us now discuss the different stages of the measurement procedure (Fig. 2.) 

a) The size-quantized levels in the dots are initially empty (lie above the chemical potentials 
in the electrodes). The dashed line shows the levels in the dots split off by the Coulomb repulsion. 

b,c) The central and the left dots are subjected to the voltage pulses with the duration r 
such that r res <r< Td yn and the central and left dots levels are tuned into the resonance and 
pulled down below the chemical potential pl in the left electrode. During the time of the order 
of T res the levels in the left and the central dots are filled by the electrons from the left electrode. 

d) Then the voltage pulse with the duration r res « r « Td yn is applied to the left dot 
pushing its level above the chemical potential pi. During the time of the order of T res the level 
in the left dot becomes empty since the electron escapes back into the left electrode. AT the 
same time the level in the central dot remains filled. On the time scale < r non one can assume 
that the electron does not remember about the electrodes and is effectively isolated, its spin 
state being determined by the left electrode state. The above described procedure results in the 
preparation at the initial moment of time (on the time scale r <C Td yn — instantly) of the ancilla 
in the state pA{t = 0). Since the density matrix for the spin-1/2 system can always be written 
in the form p = ^(-f + <r • u), the state of the electron in the central dot which tunnelled from 
the left electrode is described by the density matrix pA(t = 0) = (1/2)(I + a ■ Ul), where u^ 
is the vector describing both the direction and the magnitude of the electron spin polarization 
vector of the electrons in the left electrode. 

Then on the time scale T res <C t <C T non one can assume that the electron in the central dot 
and the spins in the gate evolve according the joint unitary dynamics 

p(t) = U(t){p A (t = 0) <g> p s (t = O))^- 1 ^), U(t) = exp(i / H int {t')dt'). 



Here p s (t = 0) is the density matrix of the quantum gate the time moment t = 0. The Hamilto- 
nian Hi n t can be easily diagonalized since it describes a finite-dimensional system (e.g., the joint 
dynamics of the electron localized in the dot together with the nucleus spin and the electron 
spin localized on it is described by the 8x8 matrix). 

The density matrix of the electron in the central dot by the times moment t after the joint 
evolution is 

PA {t) = Tr s {U(t)(p A (t = 0) <8> p s (t = OWit)- 1 } = i(J + a ■ u A (t)), 

where the vector u A (t) specifies the spin of electron localized in the central dot by the time 
moment t resulting from the interaction with the spins in the quantum gate. 

e) Detection of the electron state in the central dot is performed by measuring the current 
flowing into the right electrode. For that purpose the central and the right quantum dots are 
subjected to the voltage pulses of duration t\ (Fig. 2e) similar to those used to inject an 
electron from the left electrode to the central dot. If the time t\ is short compared with T res , the 
probability of electron escape to the right electrode is proportional to t\. Since Td yn 3> Tres, at 
the time moment t the interaction between the ancilla and the quantum gate is almost instantly 
(on the time scale characteristic of their joint dynamics) turned off. The probability of electron 
escape to the right electrode per unit time is (to within the numerical factors) 

Pr oc \T\ 2 Tr sA {p R ■ p a (t)} = \T\ 2 Tr sA {(I s ® p R )(U(t)(p s (t = 0) <g> p A (t = 0))U-\t))}, 

where p R is the electron density matrix in the right electrode, p R = -^{I + cr • Ur). Therefore, 
the probability of appearance of a current pulse in the right electrode depends on the spin state 
of the electron in the central quantum dot and is 

Pi = Cr l -\T\ 2 {l + u R -u A (t)}, (4) 

where C is a constant (the scalar product actually arises from the reduction to a single spin 
quantization axis of the two spinors describing the electron states in the central dot and the 
right electrode [20]). 

f) Finally, the voltage pulses are applied to the left and the central quantum dots whose 
magnitude and duration chosen in such a way that the electron escapes from the central dot (if 
after the previous stage there is still an electron in the central quantum dot) to the left electrode 
with unit probability. 

Let the duration of the complete cycle consisting of the stages a)-f) be tq. Then at fixed t\ 
the current flowing through the system of quantum dots will be equal to Pr • ^-. The constant 
C in Eq. (4) can be found from the current measurements for the case of parallel magnetizations 
in the left and right electrodes when the interaction with the gate is turned off (t = 0), so that 
ur • u A (t) = \u R \ • \ul\ (we assume that the magnetizations in the electrodes, |ur| and |u^|), 
are known). 

Thus, the current measurements in the outlined scheme allow one to determine the ancilla po- 
larization vector u A (t) which depends on the initial state of the gate p s before the measurement. 
Strictly speaking, finding the gate state from the measured current requires the determination of 
all three components of the vector u s appearing in the density matrix of the gate. It is obvious 
that this can only be done if the current measurements are preformed for at least three different 
combinations of the system parameters. For example, one can vary the magnetization direction 



in both electrodes and the duration of ancilla interaction with the gate. However, the problem 
of whether or not the tunnel current behavior as a function of the indicted parameters provides 
sufficient information for the complete recovery of the vector u s should be solved separately for 
each particular interaction between the quantum gate and the ancilla. 

The characteristic time of the non-resonant tunneling can be made arbitrarily large by in- 
creasing the width of the double barrier so that it imposes no restrictions. The typical time of 
the joint quantum dynamics of the gate and the electron in the central dot can be estimated 
as the typical time of the dynamics of an isolated gate which is the inverse Larmor spin pre- 
cession in the external field [8]. In the field B ~ 100 Gs (0.01 T) this time is 1/Td yn ~ 10 6 
Hz. The resonant tunneling time can well be increased up to r res ~ 10~ 9 s allowing to measure 
the current pulses on the times r < r res . To avoid the smearing of the Zeeman splitting, the 
temperature should not exceed 1 mK. Growth of the operational temperature shortens Td yn and, 
consequently, r res , resulting in the reduction of the time during which the current pulses are 
measured. Note also that both the quantum dots and silicon matrix [8] should not contain the 
isotopes with non-zero nuclear spin which prevents employment of the advanced technology of 
GaAs/GaAlAs materials and requires usage of the Si/SiGe-based systems. 
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